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Abstract 



Recent evaluations of the hadronic vacuum polarization contributions to the effective fine- 
structure constant a em (Mz) are summarized and commented. A new update based on 
00 ' corrected CMD-2 data is presented. My new estimates are Aaj&CMl) = 0.027773 ± 

0.000354 (e + e--data based) and Aa£i(M§) = 0.027664±0.000173 (via the Adler-function 
and extended use of pQCD). Prospects of further possible progress is discussed. 

1 Introduction 

Precision physics requires appropriate inclusion of higher order effects and the knowl- 
edge of very precise input parameters of the electroweak Standard Model SM. One of the 
basic input parameters is the fine structure constant which depends logarithmically on 
the energy scale. Vacuum polarization effects lead to a partial screening of the charge in 
the low energy limit (Thomson limit) while at higher energies the strength of the electro- 
magnetic interaction grows. Non-perturbative strong interaction effects (virtual hadron 
fluctuations) lead to a partially non-perturbative relationship between the very precisely 
known classical fine structure constant a and its effective value «(/z 2 ) at nonzero energy 
scales /i. Presently, the only save way to evaluate the non-perturbative hadronic contribu- 
tions is via a dispersion integral over experimental e + e~ — > hadrons data (see e.g. PUE]). 
The required relationship derives from analyticity (as a consequence of causality) and the 
optical theorem (unitarity). A drawback is that the experimental errors allow us to cal- 
culate the shift in the fine structure constant only at limited accuracy. In fact, at present, 
the corresponding uncertainty imposes one of the limiting factors in making precise SM 
predictions. While 5Aa is dominating the uncertainty S sin 2 y the experimental error 

*To appear in the Proceedings of the Mini- Workshop "ELECTROWEAK PRECISION DATA 
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of the top mass measurement 5m t is the main uncertainty contributing to 5 My/. Both 
measurements sin 2 Of and Mw yield constraints on the as yet unknown Higgs mass M#. 

The photon propagator provides a simple way to derive the concept of an effective 
charge. Including a factor e 2 and considering the renormalized ( multiplied by the appro- 
priate wave function renormalization factor Z) full photon propagator we have 

e 2 D^(q) = ff- 6 Z , + gauge terms (1.1) 

t (i + n;(g 2 )j 

which in effect means that the charge has to be replaced by a running charge 

The wave function renormalization factor Z is fixed by the condition that for q 2 — > 
one obtains the classical charge (charge renormalization in the Thomson limit). Thus the 
renormalized charge is 

^ «V) - TTgwbw) (L3) 

where, in perturbation theory, the lowest order diagram which contributes to 11^ (g 2 ) is 

nV) = 

which describes virtual creation and reabsorption of fermion pairs 

7* — ► e + e~, t + t~, uu, dd, ■ ■ — ► 7* 

in leading order. In terms of the fine structure constant a = e 2 /47r Eq. (jl.3|) reads 

J2\ _ a . a r>„ fn' t-2\ 




a{q ) = Y^Aa ; Aa = ~ Re ^ {q ) ~ U '< {0) ) ■ (L4) 

The shift Aa is large due to the large change in scale going from zero momentum to 
the Z-mass scale \x = Mz and due to the many species of fermions contributing. Zero 
momentum more precisely means the light fermion mass thresholds. 

The various contributions to the shift in the fine structure constant come from the 
leptons (lep = e, /1 and r) the 5 light quarks (u, b, s, c, and b and the corresponding 
hadrons = had) and from the top quark: 

Aa = Aa lcp + A< 5 Wd + Aa top + • • • (1.5) 

Also W^-pairs contribute at q 2 > 2M^ (see 011]). The leptonic contributions are calcu- 
lable in perturbation theory where at leading order the free lepton loops yield 



Aai e p(s) = 



£ jfc In (s/m 2 ) - § + 0(mj/s) for |s| > m 



■2 



1.6) 



0.03142 for s = Mf 



where = Jl — 4m|/s. This leading contribution is affected by small electromagnetic 
corrections only in the next to leading order. The leptonic contribution is actually known 
to three loops jSJE] at which it takes the value [M z ~ 91.19 GeV) 1 

Aai cp (Mf) ~ 314.98 x 1(T 4 . (1.7) 

In contrast, the corresponding free quark loop contribution gets substantially modified 
by low energy strong interaction effects, which cannot be obtained by perturbative QCD 
(pQCD). As already mentioned, fortunately, one can evaluate this hadronic term Aaj^ 
from hadronic e + e~- annihilation data by using a dispersion relation. The relevant once 
subtracted vacuum polarization amplitude (jl.4j) satisfies a convergent dispersion relation 
and correspondingly the shift of the fine structure constant a is given by 



where 



a(e + e — > 7* — > hadrons) . 

-A— '- — ^ = 127rImnL^(s 

cr(e+e _ — > 7* — > /j + 

Accordingly, the one particle irreducible (lpi) blob 



<had 

which is the relevant building block in our context and is given by diagrams which cannot 
be cut into two disconnected parts by cutting a single photon line, at low energies exhibits 
intermediate states like 7r°7, p, u, 4>, ■ ■ ■ , rnr, 3n, An, • ■ ■ , tttcj, titiZ, • ■ ■ , imH, ■ • • , KK, ■ ■ ■ 
(at least one hadron plus any strong, electromagnetic or weak interaction contribution) 
and the corresponding contributions are to be calculated via a dispersion relation from 
the imaginary parts which are given by the production of the corresponding intermediate 
states in e + e~-annihilation via virtual photons (at energies sufficiently below the point 
where 7 — Z interference comes into play). 

A direct evaluation of the _R(s)-data up to y/s = E cut = 5 GeV and for the T 
resonance-region between 9.6 and 13 GeV and applying perturbative QCD from 5.0 to 
9.6 GeV and for the high energy tail above 13 GeV at Mz = 91.19 GeV yields 2 : 

Aai 5 a d( M f) = 0.027773 ± 0.000354 (1.10) 
a~ l (Ml) = 128.922 ± 0.049 . 

The contributions from different energy ranges are shown in Tab. 
iFor mt - 174.3 GeV we have Aa top (M|) ~ -g^^l -6 x 1(T 5 . 

2 pQCD for calculating R(s), as worked out to high accuracy in Refs. [7]-|H], is used here only where 
it has been checked to work and converge well: in non-resonant regions at sufficiently high energies and 
sufficiently far from resonances and thresholds. I have further checked that results obtained with my own 
routines agree very well with the ones obtained via the recently published program rhad-1 .00 jl()| . 

3 Table 1 also specifies largely details of the error handling. The different energy ranges mark typical 
generation of experiments within which systematic errors are considered to be 100% correlated, while all 
errors are treated as independent for all entries of the table. 
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final state 


energy range (GeV) 


Aaf^Mf) 


(stat) (syst) 


A «Ld(- s o; 


(stat) (syst) 


XPT 


(0.28, 0.32) 


0.04 ( 


0.00) 


; o.oo) 


0.03 ( 


0.00) 


; o.oo) 


P 


(0.28, 0.81) 


26.16 ( 


0.24) 


{ 0.27) 


24.26 ( 


0.23) 


{ 0-25) 


to 


(0.42, 0.81) 


3.02 ( 


0.04) 


; 0.08) 


2.75 ( 


0.03) 


{ 0.07) 




(1.00, 1.04) 


4.74 ( 


0.07) 


I o.ii) 


4.07 ( 


0.06) 


{ 0.09) 


J/V 




11.50 ( 


0.56) 


' 0.61) 


4.06 ( 


0.19) 


{ 0.19) 


T 




1.27 ( 


0.05) 


; 0.07) 


0.07 ( 


0.00) 


; o.oo) 


hadrons 


(0.81, 1.40) 


12.92 ( 


0.13) 


' 0.52) 


11.05 ( 


0.11) 


{ 0.43) 


hadrons 


(1.40, 3.10) 


27.13 ( 


0.11) 


; o.6o) 


15.75 ( 


0.06) 


{ 0.37) 


hadrons 


(3.10, 3.60) 


5.31 ( 


0.11) 


[ o.io) 


1.90 ( 


0.04) 


[ 0.04) 


hadrons 


(3.60, 9.46) 


51.49 ( 


0.25) 


; 3.00) 


8.41 ( 


0.04) 


; 0.44) 


hadrons 


(9.46,13.00) 


18.59 ( 


0.25) 


{ 1.36) 


0.90 ( 


0.01) 


{ 0.07) 


perturb 


(13.0,oo) 


115.59 ( 


0.00) 


I 0.12) 


1.09 ( 


0.00) 


{ 0.00) 


data 
total 


(0.28,13.00) 


162.14 ( 
277.73 ( 


0.74) 
0.74) 


; 3.46) 
; 3.46) 


73.21 ( 
74.30 ( 


0.33) 
0.33) 


{ 0.80) 
{ 0.80) 



Table 1: Results for Ao;[ 5 a } d (M|) data and Aa^ d (-s ) data (^=2.5 GeV). 



Our analysis is as close to the experimental results as possible by utilizing the trape- 
zoidal rule together with PDG rules for taking weighted averages between different ex- 
periments as described in detail in [2] 

The most important ingredient of our analysis are the e + e~-data which we described 
in detail in |5j (see also jTT]) and the new data which have become available since then. 
The developments concerning the experimental data as well as some theoretical aspects 
are the following: 

• The updated results of the precise measurements of the processes e + e~ — > p — ► 7r + 7r~, 
e + e~ — > uj — ► 7T + 7r~7r° and e + e~ — >• —>■ K^Ks performed by the CMD-2 collaboration 
which have just been presented [T2] . The update appeared necessary due to an overes- 
timate of the integrated luminosity in previous analyzes 4 . The latter was published in 
2002 13J. A more progressive error estimate (improving on radiative corrections, in par- 
ticular) allowed a reduction of the systematic error from 1.4% to 0.6 % . Also some other 
CMD-2 and SND data at energies E < 1.4 GeV have become available and have been 
included. 

• Before in 2001 BES-II published their final i?-data which, in the region 2.0 GeV to 5.0 
GeV, allowed to reduce the previously huge systematic errors of about 20% to 7% [T4] . 

• After 1997 precise r-spectral functions became available [121 EH HI] which, to the 
extent that flavor SU(2){ in the light hadron sector is a symmetry, allows to obtain the 
iso-vector part of the e + e~-cross section [TB1HH]- This possibility has first been exploited 
in the present context in [11] . 

• With increasing precision of the low energy data it more and more turned out that 
we are confronted with a serious obstacle to further progress: in the region just above 
the cu-resonance, the iso-spin rotated r-data, corrected for the known iso-spin violating 

This affects in particular the leading hadronic contribution to the anomalous magnetic moment of 
the muon. 1 now obtain a)T d ^ = (695.5 ±8.6) x 10~ 10 (e + e~-data based). 
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effects, do not agree with the e + e~-data at the 10% level [20] • Before the origin of this 
discrepancy is found it will be hard to make further progress in pinning down theoretical 
uncertainties. 

• In this context iso-spin breaking effects in the relationship between the r- and the e + e~- 
data have been extensively investigated in [21]. Whatever uncertainties of the estimated 
iso-spin violations might remain, it is very unlikely that they can be made responsible for 
the observed discrepancies. 

• New results for hadronic e + e~ cross-sections are expected soon from KLOE, BABAR 
and BELLE. These experiments, running at fixed energies, are able to perform measure- 
ments via the radiative return method [221123 121] • Results presented recently by KLOE 
seem to agree very well with the final CMD-2 e + e~-data. 




Figure 1: The running of a. The "negative" E axis is chosen to indicate space-like 
momentum transfer. The vertical bars at selected points indicate the uncertainty. In 
the time-like region the resonances lead to pronounced variations of the effective charge 
(shown in the p — to and region). 

Fig. illustrates the running of the effective charge at lower energies in the space-like 
region. Typical values are Aa(5GeV) ~ 3% and Aa(Mz) ~ 6%, where about ~ 50% of 
the contribution comes from leptons and about ~ 50% from hadrons. 

An analysis similar to ours, however, using piecewise linear approximants to the non- 
resonant R(s), which are then integrated analytically, yields [23] 

Aa£d(M|) = 0.027680 ± 0.000360 (1.11) 
a-\M 2 z ) = 128.935 ± 0.049 . 

The precise choice of the lattice used for the linearization remains somewhat unclear in 
this method. 

An evaluation via the Adler-function which allows to utilize safely perturbative QCD 
for the latter at Euclidean energies above 2.5 GeV (and data at lower energies) yields (see 
below) : 

Aa£d( M f) = 0.027664 ±0.000173(0.000137] (1.12) 
ar 1 {M 2 z ) = 128.937 ± 0.024 . 
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The errors here are dominated by the QCD parameter uncertainties and are given for 
the worst case correlation [uncorrelated] case. In contrast to the so called theory-driven 
approaches, the Adler function methods is more elaborate (because it requires to transform 
data and theory to the space-like region) but allows for a dramatically better control of 
the validity of pQCD. 

There are also theory-driven approaches which utilize perturbative QCD directly for 
the evaluation of R(s) (and/or for rescaling of the normalization of the data) and assuming 
some local version of quark-hadron duality A recent result is (2H1I2Z! (see also 

Aag^Mf) = 0.027690 ±0.000180 (1.13) 
a -1 (Mf) = 128.933 ± 0.025 . 

2 ai em (s) in precision physics 

A major drawback of the partially non-perturbative relationship between a(0) and 
a(Mz) is that one has to rely on experimental data exhibiting systematic and statistical 
errors which implies a non-negligible uncertainty in our knowledge of the effective fine 
structure constant. In precision predictions of gauge boson properties this has become a 
limiting factor. Since a, and Mz are the most precisely measured parameters, they 
are used as input parameters for accurate predictions of observables like the effective weak 
mixing parameter sin 2 B/, the vector Vf and axial- vector a/ neutral current couplings, the 
W mass M w the widths Y z and Y w of the Z and the W, respectively, etc. However, for 
physics at higher energies, we have to use the effective couplings at the appropriate scale, 
for physics at the Z-resonance, for example, a(M z ) is more adequate to use than a(0). 
Of course this just means that part of the higher order corrections may be absorbed into 
an effective parameter. If we compare the precision of the basic parameters 

^ ~ 3.6 x 10- 9 ~ 1.6 -=-6.8 x 10" 4 , x 

a a(M z ) (2 I) 

^ ~ 8.6 x lO" 6 ~ 2.4 x 10~ 5 1 ' ; 

we observe that the uncertainty in a(Mz) is roughly an order of magnitude worse than 
the next best, which is the Z-mass. Future TESLA requirements are ^gg^ ~ 5.3 x 10~ 5 
(based on the assumption that a precision 5sin 2 Q f cS ~ 0.000013 (GigaZ option) and 



5M W ~ 6MeV (MegaW) may be reached) 

Let me remind the reader that Aa enters in electroweak precision physics typically 
when calculating versions of the weak mixing parameter sin 2 0j from a, and M z via 

sin 2 6; cos 2 6 i = ^— (2.2) 

V2G^M 2 Z 1- An 1 ; 

where Arj = Arrest, G^, M z ,m>H, rrif^t, m t ) includes the higher order corrections, which 
can be calculated in the SM or in alternative models. Ar has been calculated for the first 
time by A. Sirlin in 1980 (US] . In the SM, today, the Higgs boson mass m# is the only 
relevant unknown parameter and by confronting the calculated with the experimentally 
determined value of sin 2 0j one obtains important indirect constraints on the Higgs mass. 
Arj depends on the definition of sin 2 0j. The various definitions coincide at tree level and 
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hence only differ by quantum effects. From the weak gauge boson masses, the electroweak 
gauge couplings and the neutral current couplings of the charged fermions we obtain 

sin 2 e w = (2.3) 



8in2e ' = wu( l - v i)- s+v < (2 - 5) 

for the most important cases. Ari usually is written in the form 

Ar { = A« - /j(sin 2 0i) Ap + Ar ircminder (2.6) 

with a universal term Aa which affects the predictions of Mw, Afj?, A-fb-> T/, etc. The 
uncertainty 5Aa implies uncertainties 8Mw, S sin 2 Oj given by 5 

SMw 1 sm 2 Q w oi q x /\ 

~ — k— — dAa ~ 0.213 dAa (2.7) 

M w 2 cos"' Q w — sin Q w 

5sin 2 0f cos 2 6 f rA k j c* a ,„ m 

. 2n j ~ — 77^ Attp 5Aa ~ 1.54 <JAa . 2.8 

sm 2 6 / cos 2 6/ - sin" 5 6/ 

The present indirect Higgs mass "measurement" reads m H = 96tH GeV. The discrepancy 
between e + e~- and r-data based evaluations amounts to Smn ~ — 19 GeV [20] (the direct 
lower bound is itlh > 114 GeV at 95% CL while the indirect upper bound reads m# < 219 
GeV at 95% CL (1-sided). For more details we refer to [37] (in these proceedings). 

3 The r vs. e + e~ problem 

The iso-vector part of a(e + e~ — > hadrons) may be calculated by an iso-spin rotation 
from r-decay spectra, to the extend that the conserved vector current is conserved (CVC). 
The relation may be derived by comparing diagrams like: 





5 This compares with the second major source of uncertainty coming from the top mass [currently 
Smt/mt ~ 2.9 x 1CT 2 ] 

SM W cos 2 9vk Sm t Sn%t 

— — — ~ ~ 1.42b 

M w cos 2 &w - sin @w m t m t 

S sin 2 Of _ _ 2 cos 2 6/ ^£ _ _ 2 852 — 
sin 2 9 f cos 2 6/ — sin 2 0/ mt m t 
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Thus comparing t~ — > X~v T with e + e — > X° the hadronic states X~ and X° are 
approximately related by an iso-spin rotation if the states are 1 = 1 iso-vector states. 
The e + e~ cross-section is then given by 



r /=1 

r e+e- 



-v x - 



s < m T 



(3.1) 



in terms of the r spectral function vy. The r spectral function iy(s) for a given vector 
hadronic state V is defined by 



v v (s) = 



mz 



B( 



v T v~ 



1 dX 



V 



6 iKrfl 2 SW £(t 



X v ds 




(3.2) 



where |V^f| = 0.9752 ± 0.0007 j38j denotes the CKM weak mixing matrix element and 
Sew — 1-0233 ± 0.0006 accounts for electroweak radiative corrections. The spectral 
functions are obtained from the corresponding invariant mass distributions. 

Before a precise comparison is possible all kind of iso-spin breaking effects have to be 
taken into account. As mentioned earlier, this has been investigated in [21] for the most 
relevant tttc channel. Writing 



(o) 



with 



KM 



G% 



a 



\V ud \ 2 ml 



K a (s) 
KM 



dT 



T7r [ 7 ] Ryb(s) 

X 



EW 



ds 



S: 



(3.3) 



384tt 3 




~3s 



(3.4) 



the iso-spin breaking correction 



R 



IBS 



Gem(s) (3l + „ 



Ms) 



(3.5) 



include the photonic corrections (based on scalar QED), phase space corrections due to 
the — 7T° mass difference and the form-factor corrections which are dominated by the 
p — uj mixing effects. These corrections were applied in [20] and were not able to resolve 
the puzzle of the observed discrepancy (see 20] for details) 6 . 



6 The only large effect I am aware of (of order 10%) which is in the game of the comparison is a possible 
shift of the invariant mass of the pion-pairs in the p resonance region. An idea one gets if one is looking 
at the experimental p-mass values, shown in the particle data tables jSH] ("dipole shape"). If the energy 
calibration of the 7T7r-system would be to low in e + e~~ measurements or to high in r measurements by 
1% one could easily get a 10% decrease or increase in the tail, respectively. Since the p ± — p° mass 
difference as well as the difference in the widths T ± '°(p — > irir, irirj) are neither experimentally nor 
theoretically established, corresponding iso-spin violations cannot be corrected for appropriately. Note 
that the subtraction of the large and strongly energy dependent vacuum polarization effects (see Fig. ^| 
necessary for the e + e~-data, which seems to worsen the problem, is properly treated in the analysis. 
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4 Controlling pQCD via the Adler function 

In view of the increasing precision LEP experiments have achieved during the last 
few years, more accurate theoretical prediction became desirable. As elaborated in the 
introduction, one of the limiting factors is the hadronic uncertainty of Ao^ d (M|). Be- 
cause of the large uncertainties in the data, many authors advocated to extend the use of 
perturbative QCD in place of data [2E] |23|- The assumption that pQCD may be reliable 
to calculate (jl.9|) down to energies as low as 1.8 GeV seems to be supported by 

• the apparent applicability of pQCD to r physics. In fact the running of a s (m T ) — > 
a s (Mz) from the r mass up to LEP energies agrees well with the LEP value. The 
estimated uncertainty may be debated, however. 

• the smallness (see also: jXJ of non-perturbative (NP) effects if parameterized 
as prescribed by the operator product expansion (OPE) of the electromagnetic current 
correlator 

n; NP (Q 2 ) = 
+ 

+ 
+ 

where a = a s (n 2 )/n and l qfl = \n(Q 2 /fi 2 ). < —GG > and < m q qq > are the scale- 
invariantly defined condensates. 

Progress in pQCD here comes mainly from |4L)| . In addition an exact two-loop calcula- 
tion of the renormalization group (RG) in the background field MOM scheme (BF-MOM) 
is available [33] . This allows us to treat "threshold effects" closer to physics than in the MS 
scheme. The BF-MOM scheme respects the QCD Slavnov- Taylor identities (non-Abelian 
gauge symmetry) but in spite of that is gauge parameter (£) dependent 7 . 

In Ref. |32] a different approach of pQCD improvement was proposed, which relies on 
the fact that the vacuum polarization amplitude II (q 2 ) is an analytic function in q 2 with 
a cut in the s-channel q 2 = s > at s > 4m 2 and a smooth behavior in the i-channel 
(space-like or Euclidean region). Thus, instead of trying to calculate the complicated 
function R(s), which obviously exhibits non-perturbative features like resonances, one 
considers the simpler Adler function in the Euclidean region. In |32] the Adler function 
was investigated and pQCD was found to work very well above 2.5 GeV, provided the 
exact three-loop mass dependence was used (in conjunction with the background field 
MOM scheme). The Adler function may be defined as a derivative 

D{- 8 ) = -(12vr 2 ) s = ^4 Aahad(s) (4.2) 
as a as 

In applications considered below all numerical calculations have been performed in the "Landau 

gauge" £ = 0. 



Ana 



q=u 



E Q 2 q N ( 
ii 



1 

12 



11 \ < ^GG > 



-a 



a 

1 + 3 + 



_ - h 

2 ^ w ' a 



18 "J Q 4 
o ^ < m q qq > 



4 M 257 1 

— a + -Cq l aLl 

27 V3 S 486 3 m 



a 



Q 4 

q'=u,d 



< m q ,q'q' > 

Q 1 



(4.1) 
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of (|1.9|) which is the hadronic contribution to the shift of the fine structure constant. It 
is represented by 

in terms of the experimental e + e~-data. The standard evaluation (J2j) of (|4.H|) then yields 
the non-perturbative "experimental" Adler function, as displayed in Fig. El for the lower 
energies where it becomes non-perturbative. 



4.0: 


/ / / / 
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2-5-j 
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X 
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1.0: 




pQCD 
pQCD 
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n f =5 3-loop 


n 
\\ 
\\ 


0.5 : 


■-M-- 


pQCD 


n f =5 3-loop * NP 


\\ 
li 


o.O: 




data 


ncl. BESII 01/02, CMD-2 03 


\ 



-14.0 -12.0 -10.0 -8.0 -6.0 -4.0 -2.0 0.0 

Q (GeV) 



Figure 2: Adler function: theory vs. experiment |4*2"] . 

For the pQCD evaluation it is mandatory to utilize the calculations with massive 
quarks which are available up to three-loops jlU]- The four- loop corrections are known in 
the approximation of massless quarks [7] . The outcome of this analysis is pretty surprising 
and is shown in Fig- El For a discussion we refer to the original paper |4"2*j . The result was 
obtained using the background-field MOM renormalization scheme, mentioned before. 
In the transition from the MS to the MOM scheme we adapt the rescaling procedure 
described in jHj, such that for large \x 

a s ((x n) 2 ) = a s (/i 2 ) + + 0(a 3 s ) . 

This means that xq is chosen such that the couplings coincide to leading and next -to- 
leading order at asymptotically large scales. Numerically we find xq — 2.0144. Due to this 
normalization by rescaling the coefficients of the Adler-function remain the same in both 
schemes up to three-loops. In the MOM scheme we automatically have the correct mass 
dependence of full QCD, i.e., we have automatic decoupling and do not need decoupling 
by hand and matching conditions like in the MS scheme. For the numerical evaluation 
we use the pole quark masses |HB] m c = 1.55GeV, m& = 4.70GeV, m t = 173.80GeV and 
the strong interaction coupling a^^(M z ) = 0.120 ± 0.003. For further details we refer 
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According to (|4.2|) . we may compute the hadronic vacuum polarization contribution 
to the shift in the fine structure constant by integrating the Adler function. In the region 
where pQCD works fine we integrate the pQCD prediction, in place of the data. We thus 
calculate in the Euclidean region 



A^(-MI) = [A^(-MI) - Aa&(-* )] pQOD + Aa&C-so)** . (4.4) 
A save choice is so = (2.5 GeV) 2 where we obtain 

A "had(-so) data = 0.007430 ± 0.000087 (4.5) 

from the evaluation of the dispersion integral (|1.9|) . With the results presented above we 
find 

Aa£Jd(-Aff ) = 0.027626 ± 0.000087 ± 0.000149[0.000101] (4.6) 

for the Euclidean (t-channel) effective fine structure constant. The second error comes 
from the variation of the pQCD parameters. In square brackets the error if we assume 
the uncertainties from different parameters to be uncorrelated. The uncertainties coming 
from individual parameters are listed in the following table (masses are the pole masses): 



parameter range pQCD uncertainty total error 



a s 0.117... 


0.123 


0.000051 


0.000155 


m c 1.550 ... 


1.750 


0.000087 


0.000170 


m b 4.600 ... 


4.800 


0.000011 


0.000146 


m t 170.0 ... 


180.0 


0.000000 


0.000146 


all correlated 




0.000149 


0.000209 


all uncorrelated 




0.000101 


0.000178 



The largest uncertainty is due to the poor knowledge of the charm mass. I have taken 
errors to be 100% correlated. The uncorrelated error is also given in the table. 
Comments: 

• Contributions to the Adler function up to three-loops all have the same sign and are 
substantial. Four- and higher-orders could still add up to non-negligible contribution. 
An error for missing higher order terms is not included. The scheme dependence MS 
versus background field MOM has been discussed in Ref. |4*Tj . 

• The effective fine structure constant in the time-like region (s-channel), as required for 
e + e~-collider physics, may be obtained from the Euclidean one by adding the difference 

A = Aa^(Ml) - Aag^-Mf) = 0.000038 ± 0.000005 , (4.7) 

which may be calculated perturbatively or directly from the "non-perturbative" 8 disper- 
sion integral. It accounts for the i7r-terms 

M-V/V) = H\q 2 /A) + i* 

8 Since we utilize pQCD for the high energy tail in the dispersion integral, A(s) for large s is dominated 
by the tail and thus in fact is perturbative. 
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from the logs. 

• One may ask the question whether these terms should be resummed at all, i.e., included 
in the running coupling. Usually such terms tend to cancel against constant rational terms 
which are not included in the renormalization group (RG) evolution. It should be stressed 
that the Dyson summation (propagator bubble summation) in general is not a systematic 
resummation of leading, sub-leading etc. terms as the RG resummation is. 

It is worthwhile to stress here that the running coupling is not a true function of q 2 
(or even an analytic function of q 2 ) but a function of the RG scale fi 2 . The coupling as 
it appears in the Lagrangian in any case must be a constant, albeit a /^-dependent one, 
if we do not want to end up in conflict with basic principles of quantum field theory. 
The effective identification of /i 2 with a particular value of q 2 must be understood as a 
subtraction (reference) point. 




12.GeV -oo to -2.5 GeV 



Figure 3: Comparison of the distribution of contributions and errors (shaded areas scaled 
up by 10) in the standard (left) and the Adler function based approach (right), respec- 
tively. 

Since A Eq. (I4.7J1 is small we may include it in the resummation without further 
worrying and thus obtain 

Aa£2i(Aff) = 0.027664 ± 0.000173[0.000137] . (4.8) 

The alternative evaluation by the Euclidean approach is compared with the standard 
evaluation in Tab.[TJ The two methods (standard vs. Euclidean) of evaluating Aa^ d (M|) 
are also compared in Fig. |3] 

Our alternative procedure to evaluate Aa^ d (— M§) in the Euclidean region has several 
advantages as compared to other approaches used so far: The virtues of our analysis are 
the following: 

• no problems with the physical threshold and resonances 
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• pQCD is used only in the Euclidean region and not below 2.5 GeV. For lower scales 
pQCD ceases to describe properly the functional dependence of the Adler function 02] 
(although the pQCD answer remains within error bands down to about 1.6 GeV). 

• no manipulation of data must be applied and we need not refer to global or even local 
duality. That power corrections of the type Eq. (|4.1|) are negligible has been known for 
a long time. This, however, does not proof the absence of other kind of non-perturbative 
effects. Therefore our conservative choice of the minimum Euclidean energy seems to be 
necessary. 

• According to Tab.^our non-perturbative "remainder" Aq^(- So) is mainly sensitive 
to low energy data, which changes the chances of possible future experimental improve- 
ment dramatically, as illustrated in Fig. El 

While the uncertainties to Aa|^ d (M|) in the standard approach are coming essentially 
from everywhere below My , which would make a new scan over all energies for a precision 
measurement of Ohad = a(e + e~ — > 7* — ► hadrons) unavoidable, the new approach leads 
to a very different situation. The uncertainty of Aa^j(-s ) is completely dominated by 
the uncertainties of data below Mj/^ and thus new data on o"h a d are only needed below 
about 3.6 GeV which could be covered by a tunable "r-charm facility". 

5 Status and outlook 

Recent result obtained by different authors for the hadronic contributions to a em (M z ) 
are in fairly good agreement. The estimated uncertainties vary substantially, depending 
mainly on additional theoretical assumptions made in the analyzes. Tabled compares our 
results with results obtained by other authors which obtain smaller errors because they 
are using pQCD in a less controlled manner. 





) SAa 


5 sin 2 6/ 


5M W 


Method 


Rcf. 


0.02800 


0.00065 


0.000232 


11.1 


data < 12. GeV 


13 


0.027773 


0.000354 


0.000126 


6.1 


[2] + new data CMD & BES 


flTEH) 


0.027664 


0.000173 


0.000062 


3.0 


J Euclidean > 2.5 GeV 


(|T~T2l 


0.027680 


0.000360 


0.000128 


6.2 


data < 12. GeV 




0.02777 


0.00017 


0.000061 


2.9 


data < 1.8 GeV 


EDI 


0.02763 


0.00016 


0.000057 


2.7 


data < 1.8 GeV 


|32] 


0.027690 


0.000180 


0.000064 


3.1 


scaled data, pQCD 2.8-3.7, 5-oo 


I2E1I2ZI 




0.00007 


0.000025 


1.2 


5 a < 1% up to J /if) 






0.00005 


0.000018 


0.9 


5a < 1% up to T 




world 


average 


0.000170 


23.0 


LEPEWWG 2003 





Table 2: AajJ^Mf) and its uncertainties in different evaluations and their contribution 
to the errors of sin 2 Of and My^ according to (|2.8|) . Two entries show what can be reached 
by increasing the precision of cross section measurements to 1%. SMyy in MeV. 



The reduction of theoretical and experimental uncertainties must go on to cope with 
the increased energy and luminosity at future colliders (like TESLA). Ideally a reduction 
of the errors in a em (Mz) by about a factor 5 should be achieved. Such progress is equally 
important for future precision experiments at lower energies. 
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In order to be able to rely more on pQCD the required reduction of errors can be 
achieved only if the precision of QCD parameters improves accordingly. On the one 
hand this proceeds along the traditional perturbative QCD vs. experimental data line 
(see |H] and references therein), on the other hand, lattice QCD calculations will become 
of increasing importance, in this context. A lot has been achieved in this direction already 
in recent years 

However, equally important, experimental efforts must go on in measuring a(e + e~ — > 
hadrons) at the 1% level up to energies 3.6 GeV. As most of the existing facilities which 
are able to measure cr(e + e~ — > hadrons) have approved or discuss upgrade programs, we 
are confident that the progress needed actually will take place. 
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